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In these proceedings a recent calculation of the last missing piece of the two- loop 0(aa s ) 
C^l ' corrections to yti vertex at the ti threshold due to the exchange of a W boson and 

, a gluon is summarised. The calculation constitutes a building block of the top quark 

■ threshold production cross section at electron positron colliders. 

00 ' 1 Introduction 

The measurement of the total ti cross section constitutes a major part of the top physics 
f-H \ program at a future e + e~ International Linear Collider (ILC). Due to the large top quark 

mass and the large width T t 3> ^qcd, non perturbative effects are strongly suppressed at 
energies around the threshold. Therefore it is possible to determine top quark properties 
like the mass mt and the width T t , but also the strong coupling a s with high precision. In 
C"| \ Rcf. P] it has been shown that an uncertainty of below 100 MeV can be obtained for m t 

from a threshold scan of the cross section. 

The feasibility of such high-precision measurements requires a theory prediction of the 
J> ■ total cross section a(e + e~ — ► ti) with high accuracy (preferably 8a /a < 3%). The study of 

f — | ti threshold production is performed in the framework of non-relativistic QCD (NRQCD) [3] 

which separates the hard and soft scales involved in the process. The next-to-next-to-leading 
order (NNLO) calculation [JJ turned out to be as large as the NLO one. Only the (partial) 
NNNLO [6] results show a good convergence of the perturbative expansion and a strong 
reduction of the scale dependence from NNLO to NNNLO can be observed. The remaining 
uncertainty turns out to be of the order of 10%. A similar value is obtained from the 
| approach based on the resummation of logarithmically enhanced terms which have been 

considered in Refs. [7J [5]. 

In order to reach a theory goal of 8a j a < 3% it is necessary to include in the prediction 
next to the one-loop electroweak corrections, which are known since quite some time [9] (see 
also (TO] ) , also higher order effects. The evaluation of 0(aa s ) corrections has been started 
in Ref. [TT], where the two-loop mixed electroweak and QCD corrections to the matching 
coefficient of the vector current has been computed due to a Higgs or Z boson exchange 
in addition to a gluon. In these proceedings the calculation of 0(aa s ) for the two-loop 
vertex diagrams mediated by a W boson and gluon exchange is summarised [12]. This 
result completes the vertex corrections of order aa s — a building block for the top quark 
production cross section. Assuming the (numerically well justified) power counting a ~ a 2 s 
one can see that these corrections are formally of NNNLO. 

In order to complete the matching corrections of order aa s also the two-loop box dia- 
grams contributing to e + e~ — > ti have to be considered. Actually, only the proper com- 
bination of the box, vertex and self-energy contributions (the latter can, e.g., be found in 
Refs. [Ml [14]) forms a gauge independent set. 
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2 Threshold cross section 



Within the framework of NRQCD [3J the total cross section for the top quark production 
can be cast in the form 

R(e+e R ^tiX) = ^-lm[(h R , v ) 2 H v + (h R , A ) 2 H A ], (1) 

where s is the square of the centre-of-mass energy and i?(e^e^ — ► tiX) is the cross section 
normalised to cr(e + e~ — > /i+jii - ) = (47ra 2 )/(3s). For illustration we consider in Eq. |T]) 
left-landed positrons and right-handed electrons. For ete^ in the initial state a similar 
expression is obtained by replacing R by L in Eq. ([I]) . Note that in the SM the initial states 
e^e^ and eje^ are suppressed by a factor (m e /Mw) 2 ~ 10~ 10 and are thus negligible. 
We denote hay and Jir^a by hclicity amplitudes which absorb the matching coefficients 
representing the effective coupling of the effective operators. They take care of the hard 
part of the reaction. The first subscript of h refers to helicity of the electron, and the second 
one to the vector (Jy = xp^^ip) or axial-vector coupling (J^ = ^'7^75^) of the gauge bosons 
to the top quark current. The bound-state dynamics is contained in the so-called hadronic 
part which is denoted by Hy and Ha in Eq. (fTjl. 

In the cross section formula (p} "Im" refers to those cuts which correspond to the tiX 
final state. This means that we have to select special cuts which correspond to the final 
state we are interested in. This requires a dedicated study incorporating the experimental 
setup. For the one-loop electroweak correction this treatment was performed in |15j . In [12j 
we did not pursue this problem further. 

The evaluation of Hy and Ha requires to integrate out the low-energy modes of QCD, 
the soft, potential and ultrasoft gluons contained in NRQCD [TCI I17j. For the top quark 
system this can be done perturbatively. In a first step one integrates out the soft and 
potential gluons which results in the effective field theory Potential NRQCD [T8l [19] . The 
corresponding Lagrangian is known to 
NNNLO [5U|F1 To integrate out nonrela- 
tivistic top and anti-top quark fields the 
Rayleigh-Schrodinger perturbation theory 
can be applied as was initiated in Ref. [22] 
and performed to NNNLO for Hy and Ha 
in Refs. [5] and [53], respectively. Inte- 
grating out the ultrasoft gluon was com- 
pleted recently in Ref. [5]. For the details 
of these steps we refer to the original pa- 
pers and references cited therein (see also 
Rcfs. |21 [33 [2SI ED |2g|). 

In these proceedings we report on a 
recent publication, where corrections of 
0{aa s ) to hi p were computed. Some sam- 
ple vertex diagrams which have to be eval- 
uated at threshold are shown in Fig. [TJ Let 
us mention that the corresponding two-loop 

a The only missing constant in Ref. 1201 is related to the three-loop static potential where recently the 
fcrmion corrections became available |21j . 
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Figure 1: Sample Feynman diagrams con- 
tributing to the matching coefficient of the 
vector current at order a (a), a s (b) and aa s 
(c)-(f). The thick (thin) straight lines repre- 
sent top (bottom) quarks, wavy lines stand for 
W bosons and the curly ones for gluons. 
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QCD corrections have been presented in Refs. [551 [3D] an d the three-loop corrections induced 
by a light quark loop in Ref . [21] ■ 



3 Calculation of the vertex diagrams 



The two-loop diagrams which have to be evaluated are generated with QGRAF [32] and further 
processed with q2e and exp [3"31 [51] . The reduction of the integrals is performed with the 
program crusher [35 which implements the Laporta algorithm [361 137] . We arrive at 
29 master integrals (MI). Some Mis factorise into one-loop integrals or contain only one 
dimcnsionful scale. Most of these integrals are available in the literature and can, e.g., be 
found in Refs. [Ml HE SOI EH [TQ [42] . 

The two-scale Mis are evaluated via an expansion in z = M^/m 2 . A promising method 
to obtain the results is based on differential equations (see Ref. [43] for a recent review) 
which provide the expansion in an automatic way once the initial conditions are specified. 
With the help of the ansatz 



MI 



(2) 



the differential equations can be expanded in e and z. As a result they reduce to algebraic 
equations for the coefficients Cijk, which can be solved trivially. In every order in e there is 
one constant Cijk which can not be determined with this procedure. It is obtained from the 
initial condition at z = 0. Unfortunately, we could not get analytic results for five coefficients 
in the e-expansion of four Mis at z = 0. We calculated these coefficients using the Mellin- 
Barnes method (see, e.g., Ref [13]) where we used the program packages AMBRE [35] and 

mb jug. 

The Mellin-Barnes representation for a given integral 
is not unique. In particular it might happen that the 
convergence of the resulting numerical integration turns 
out to be good in one case whereas a poor convergence 
is observed in other cases. The crucial quantity in this 
respect is the asymptotic behaviour of the T function for 
large imaginary part which is given by 



T(a ± ib) 



2ne 
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where the first two exponential factors lead to oscillations. 
Let us discuss this in more detail for the Mellin-Barnes 
representation of one of the Mis, which is depicted in 
Fig. El 



Ml = 



Figure 2: Solid (dotted) lines 
correspond to propagators with 
mass m t (0). Single external 
lines are on the mass shell with 
mass m t whereas double exter- 
nal line have mass 2m t . 
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^ r(l - e)r(-e - zi + l)r(-zi) [T(-2e - z 2 + l)] 2 
X r(-3e -zi + 2)r(-2e -z x + 2)I\-4e - 2z 2 + 2) 
xr(-4e - Zl -z 2 + 2)T(z! - z 2 )r(e + z 2 )r(2e + z 2 ) , 



(4) 



where je = 0.577216 . . . and the contour of integration is chosen in such a way that the poles 
of the r functions with +Zi are separated from the poles of the T functions with — Zj. Using 
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the package MB we can expand the integrand in e. For the finite contribution this leads to a 
sum of an analytic part, a one-dimensional Mcllin-Barncs integral and a two-dimensional one. 
The latter correspond to the integral in Eq. (f4| for e = 0. If we insert in this expression the 
asymptotic behaviour for the T functions as given in Eq. (|3|) one observes that the integrand 
of the two-dimensional integral falls off exponentially, except for Im(z2) = 0,lm(zi) < 00 
On this line the drop-off only shows a power-law behaviour which is dictated by the last 
factor of Eq. ([3]) . In our particular case the drop-off turns out to be extremely slow for the 
integration contour chosen by MB which corresponds to Re(zi) = —1/4 and Re(z2) = — 1/2. 
Thus it is hard to get an accurate result by the numerical integration since a highly oscillating 
functions has to be integrated. A closer look to the fall-off behaviour in Eq. shows that 
it is possible to improve the drop-off for Im^i) — » — oo by taking residues of the integrand 
in zi and thus shifting the integration contour for more and more to positive values for 
Re(z2). In this way the integrand becomes well-behaved and can be integrated numerically 
with sufficiently high precision. 

It has already been pointed out in [46j that for certain kinematical configurations the 
Mcllin-Barnes integrals exhibit poor convergence behaviour. We have shown that at least 
for the threshold integrals needed in our calculation, it is possible to choose the integration 
paths in a way to make the numerical integration possible. One may hope that with the 
approach described here it will turn out to be possible to use Mcllin-Barncs integration for 
other troublesome integrals, too. 

The other integrals which were solved using the Mellin-Barnes method show similar 
properties than the MI above. In all cases it is possible to end up with integrals which could 
be integrated numerically. The accuracy for the finite part of these integrals is sufficient to 
obtain the final result with four significant digits. 

Note that contrary to the default settings of MB we do not use Vegas for the multidi- 
mensional numerical integrations. Instead we use Divonne which is available from the Cuba 
library [47] . For the integrals we have considered it leads to more accurate results using less 
CPU time. 

We have performed an independent check of the initial conditions for all the Mis employ- 
ing the method of sector decomposition. In particular we used the program FIESTA [48j . 

4 Results 

The hard part of the cross section in (fTJ) at tree level is given by 

hTv/A = Q e Q t + S -^^ with = (5) 

where the /?/ is the coupling of a fermion (/ = e, t) with helicity / to the Z boson [T^]. We 
denote by T A the contribution of the sum of all onc-particle-irrcduciblc diagrams to the jti 
vertex and parametrise the radiative corrections in the form 

fi = Q t + f^ ^+f^ 1 '°)+f^ 1 ' 1) , (6) 

where the hat denotes renormalised quantities and the subscripts indicate corrections 
of 0(a l a{). Substituting for the Qt in Eq. © and retaining the relevant orders in 

b Note that for these values of z\ the exponential factor in the integrand of {4} increases exponentially. 
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the electroweak and strong couplings leads to the corrections to the helicity amplitudes, 
h^y.h^y 1 and h^'y . We further decompose t\ (and similarly the quantities on the 
right-hand side of Eq. ©) according to the contributions from the Higgs, Z and W boson 
exchanges: 

fyt = T*A,H + ^ A,Z + ^A,W ■ (?) 
The results for the electroweak corrections read EJ [H] 

21.1 x 10~ 3 (10.6 x 1CT 3 ) for M H = 120 (200) GeV , 
3.0 x 10 -3 , 
1.7 x 10~ 3 , 

-17.6 x 10~ 3 (-6.6 x 10~ 3 ) for M H = 120 (200) GeV , 
0.2 x 10~ 3 , 

-1.0 xlO" 3 . (8) 

One observes quite small corrections from the W and Z boson induced contributions. From 
Eq. ([8]) one can read off that relatively big one-loop effects are obtained for light Higgs 
boson masses. However, there is a strong cancellation between the one- and two-loop terms 
resulting in corrections which have the same size as the sum of the one- and two-loop 
contributions of the W and Z boson diagrams. In general moderate effects are observed 
suggesting that in the electroweak sector perturbation theory works well, which is in contrast 
to the pure QCD corrections. Let us in the end mention that Jy 1 leads to a correction 
of 0.1% to R. 
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